Synovial fibroblast gene expression is associated with sensory nerve growth and pain in rheumatoid arthritis.

It has been presumed that rheumatoid arthritis (RA) joint pain is related to inflammation in the synovium; however, recent studies reveal that pain scores in patients do not correlate with synovial inflammation. We developed a machine-learning approach (graph-based gene expression module identification or GbGMI) to identify an 815-gene expression module associated with pain in synovial biopsy samples from patients with established RA who had limited synovial inflammation at arthroplasty. We then validated this finding in an independent cohort of synovial biopsy samples from patients who had early untreated RA with little inflammation. Single-cell RNA sequencing analyses indicated that most of these 815 genes were most robustly expressed by lining layer synovial fibroblasts. Receptor-ligand interaction analysis predicted cross-talk between human lining layer fibroblasts and human dorsal root ganglion neurons expressing calcitonin gene-related peptide (CGRP+). Both RA synovial fibroblast culture supernatant and netrin-4, which is abundantly expressed by lining fibroblasts and was within the GbGMI-identified pain-associated gene module, increased the branching of pain-sensitive murine CGRP+ dorsal root ganglion neurons in vitro. Imaging of solvent-cleared synovial tissue with little inflammation from humans with RA revealed CGRP+ pain-sensing neurons encasing blood vessels growing into synovial hypertrophic papilla. Together, these findings support a model whereby synovial lining fibroblasts express genes associated with pain that enhance the growth of pain-sensing neurons into regions of synovial hypertrophy in RA.

Deconstruction of rheumatoid arthritis synovium defines inflammatory subtypes.

Rheumatoid arthritis is a prototypical autoimmune disease that causes joint inflammation and destruction1. There is currently no cure for rheumatoid arthritis, and the effectiveness of treatments varies across patients, suggesting an undefined pathogenic diversity1,2. Here, to deconstruct the cell states and pathways that characterize this pathogenic heterogeneity, we profiled the full spectrum of cells in inflamed synovium from patients with rheumatoid arthritis. We used multi-modal single-cell RNA-sequencing and surface protein data coupled with histology of synovial tissue from 79 donors to build single-cell atlas of rheumatoid arthritis synovial tissue that includes more than 314,000 cells. We stratified tissues into six groups, referred to as cell-type abundance phenotypes (CTAPs), each characterized by selectively enriched cell states. These CTAPs demonstrate the diversity of synovial inflammation in rheumatoid arthritis, ranging from samples enriched for T and B cells to those largely lacking lymphocytes. Disease-relevant cell states, cytokines, risk genes, histology and serology metrics are associated with particular CTAPs. CTAPs are dynamic and can predict treatment response, highlighting the clinical utility of classifying rheumatoid arthritis synovial phenotypes. This comprehensive atlas and molecular, tissue-based stratification of rheumatoid arthritis synovial tissue reveal new insights into rheumatoid arthritis pathology and heterogeneity that could inform novel targeted treatments.

Machine Learning Reveals Synovial Fibroblast Genes Associated with Pain Affect Sensory Nerve Growth in Rheumatoid Arthritis.

It has been presumed that rheumatoid arthritis (RA) joint pain is related to inflammation in the synovium; however, recent studies reveal that pain scores in patients do not correlate with synovial inflammation. We identified a module of 815 genes associated with pain, using a novel machine learning approach, Graph-based Gene expression Module Identification (GbGMI), in samples from patients with longstanding RA, but limited synovial inflammation at arthroplasty, and validated this finding in an independent cohort of synovial biopsy samples from early, untreated RA patients. Single-cell RNA-seq analyses indicated these genes were most robustly expressed by lining layer fibroblasts and receptor-ligand interaction analysis predicted robust lining layer fibroblast crosstalk with pain sensitive CGRP+ dorsal root ganglion sensory neurons. Netrin-4, which is abundantly expressed by lining fibroblasts and associated with pain, significantly increased the branching of pain-sensitive CGRP+ neurons in vitro . We conclude GbGMI is a useful method for identifying a module of genes that associate with a clinical feature of interest. Using this approach, we find that Netrin-4 is produced by synovial fibroblasts in the absence of inflammation and can enhance the outgrowth of CGRP+ pain sensitive nerve fibers. One Sentence Summary: Machine Learning reveals synovial fibroblast genes related to pain affect sensory nerve growth in Rheumatoid Arthritis addresses unmet clinical need.

Computer Vision Analysis of Rheumatoid Arthritis Synovium Reveals Lymphocytic Inflammation Is Associated With Immunoglobulin Skewing in Blood.

OBJECTIVE: We sought to develop computer vision methods to quantify aggregates of cells in synovial tissue and compare these with clinical and gene expression parameters. METHODS: We assembled a computer vision pipeline to quantify five features encompassing synovial cell density and aggregates and compared these with pathologist scores, disease classification, autoantibody status, and RNA expression in a cohort of 156 patients with rheumatoid arthritis (RA) and 149 patients with osteoarthritis (OA). RESULTS: All five features were associated with pathologist scores of synovial lymphocytic inflammation (P < 0.0001). Three features that related to the cells per unit of tissue were significantly increased in patients with both seronegative and seropositive RA compared with those with OA; on the other hand, aggregate features (number and diameter) were significantly increased in seropositive, but not seronegative, RA compared with OA. Aggregate diameter was associated with the gene expression of immunoglobulin heavy-chain genes in the synovial tissue. Compared with blood, synovial immunoglobulin isotypes were skewed from IGHM and IGHD to IGHG3 and IGHG1. Further, patients with RA with high levels of lymphocytic infiltrates in the synovium demonstrated parallel skewing in their blood with a relative decrease in IGHGM (P < 0.002) and IGHD (P < 0.03) and an increase in class-switched immunoglobulin genes IGHG3 (P < 0.03) and IGHG1 (P < 0.002). CONCLUSION: High-resolution automated identification and quantification of synovial immune cell aggregates uncovered skewing in the synovium from naïve IGHD and IGHM to memory IGHG3 and IGHG1 and revealed that this process is reflected in the blood of patients with high inflammatory synovium.

Machine learning identification of thresholds to discriminate osteoarthritis and rheumatoid arthritis synovial inflammation.

BACKGROUND: We sought to identify features that distinguish osteoarthritis (OA) and rheumatoid arthritis (RA) hematoxylin and eosin (H&E)-stained synovial tissue samples. METHODS: We compared fourteen pathologist-scored histology features and computer vision-quantified cell density (147 OA and 60 RA patients) in H&E-stained synovial tissue samples from total knee replacement (TKR) explants. A random forest model was trained using disease state (OA vs RA) as a classifier and histology features and/or computer vision-quantified cell density as inputs. RESULTS: Synovium from OA patients had increased mast cells and fibrosis (p < 0.001), while synovium from RA patients exhibited increased lymphocytic inflammation, lining hyperplasia, neutrophils, detritus, plasma cells, binucleate plasma cells, sub-lining giant cells, fibrin (all p < 0.001), Russell bodies (p = 0.019), and synovial lining giant cells (p = 0.003). Fourteen pathologist-scored features allowed for discrimination between OA and RA, producing a micro-averaged area under the receiver operating curve (micro-AUC) of 0.85±0.06. This discriminatory ability was comparable to that of computer vision cell density alone (micro-AUC = 0.87±0.04). Combining the pathologist scores with the cell density metric improved the discriminatory power of the model (micro-AUC = 0.92±0.06). The optimal cell density threshold to distinguish OA from RA synovium was 3400 cells/mm2, which yielded a sensitivity of 0.82 and specificity of 0.82. CONCLUSIONS: H&E-stained images of TKR explant synovium can be correctly classified as OA or RA in 82% of samples. Cell density greater than 3400 cells/mm2 and the presence of mast cells and fibrosis are the most important features for making this distinction.

Oral mucosal breaks trigger anti-citrullinated bacterial and human protein antibody responses in rheumatoid arthritis.

Periodontal disease is more common in individuals with rheumatoid arthritis (RA) who have detectable anti-citrullinated protein antibodies (ACPAs), implicating oral mucosal inflammation in RA pathogenesis. Here, we performed paired analysis of human and bacterial transcriptomics in longitudinal blood samples from RA patients. We found that patients with RA and periodontal disease experienced repeated oral bacteremias associated with transcriptional signatures of ISG15+HLADRhi and CD48highS100A2pos monocytes, recently identified in inflamed RA synovia and blood of those with RA flares. The oral bacteria observed transiently in blood were broadly citrullinated in the mouth, and their in situ citrullinated epitopes were targeted by extensively somatically hypermutated ACPAs encoded by RA blood plasmablasts. Together, these results suggest that (i) periodontal disease results in repeated breaches of the oral mucosa that release citrullinated oral bacteria into circulation, which (ii) activate inflammatory monocyte subsets that are observed in inflamed RA synovia and blood of RA patients with flares and (iii) activate ACPA B cells, thereby promoting affinity maturation and epitope spreading to citrullinated human antigens.

Factors Associated With Increasing Length of Stay for Rheumatoid Arthritis Patients Undergoing Total Hip Arthroplasty and Total Knee Arthroplasty.

Background: Total hip arthroplasty (THA) and total knee arthroplasty (TKA) are cost-effective procedures that decrease pain and improve health-related quality of life for patients with advanced symptomatic arthritis, including rheumatoid arthritis (RA). Patients with RA have a longer length of stay (LOS) after THA or TKA than patients with osteoarthritis, yet the factors contributing to LOS have not been investigated. Purpose: We sought to identify the factors contributing to LOS for patients with RA undergoing THA and TKA at a single tertiary care orthopedic specialty hospital. Methods: We retrospectively reviewed data from a prospectively collected cohort of 252 RA patients undergoing either THA or TKA. Demographics, RA characteristics, medications, serologies, and disease activity were collected preoperatively. Linear regression was performed to explore the relationship between LOS (log-transformed) and possible predictors. A multivariate model was constructed through backward selection using significant predictors from a univariate analysis. Results: Of the 252 patients with RA, 83% were women; they had a median disease duration of 14 years and moderate disease activity at the time of arthroplasty. We had LOS data on 240 (95%) of the cases. The mean LOS was 3.4 ± 1.5 days. The multivariate analysis revealed a longer LOS for RA patients who underwent TKA versus THA, were women versus men, required a blood transfusion, and took preoperative opioids. Conclusion: Our retrospective study found that increased postoperative LOS in RA patients undergoing THA or TKA was associated with factors both non-modifiable (type of surgery, sex) and modifiable (postoperative blood transfusion, preoperative opioid use). These findings suggest that preoperative optimization of the patient with RA might focus on improving anemia and reducing opioid use in efforts to shorten LOS. More rigorous study is warranted.

Dynamic spectrum of ectopic lymphoid B cell activation and hypermutation in the RA synovium characterized by NR4A nuclear receptor expression.

Ectopic lymphoid structures (ELS) can develop in rheumatoid arthritis (RA) synovial tissue, but the precise pathways of B cell activation and selection are not well understood. Here, we identify a synovial B cell population characterized by co-expression of a family of orphan nuclear receptors (NR4A1-3), which is highly enriched in RA synovial tissue. A transcriptomic profile of NR4A synovial B cells significantly overlaps with germinal center light zone B cells and an accrual of somatic hypermutation that correlates with loss of naive B cell state. NR4A B cells co-express lymphotoxins α and ß and IL-6, supporting functions in ELS promotion. Expanded and shared clones between synovial NR4A B cells and plasma cells and the rapid upregulation with BCR stimulation point to in situ differentiation. Together, we identify a dynamic progression of B cell activation in RA synovial ELS, with NR4A transcription factors having an important role in local adaptive immune responses.

Rheumatoid Arthritis Synovial Inflammation Quantification Using Computer Vision.

OBJECTIVE: We quantified inflammatory burden in rheumatoid arthritis (RA) synovial tissue by using computer vision to automate the process of counting individual nuclei in hematoxylin and eosin images. METHODS: We adapted and applied computer vision algorithms to quantify nuclei density (count of nuclei per unit area of tissue) on synovial tissue from arthroplasty samples. A pathologist validated algorithm results by labeling nuclei in synovial images that were mislabeled or missed by the algorithm. Nuclei density was compared with other measures of RA inflammation such as semiquantitative histology scores, gene-expression data, and clinical measures of disease activity. RESULTS: The algorithm detected a median of 112,657 (range 8,160-821,717) nuclei per synovial sample. Based on pathologist-validated results, the sensitivity and specificity of the algorithm was 97% and 100%, respectively. The mean nuclei density calculated by the algorithm was significantly higher (P < 0.05) in synovium with increased histology scores for lymphocytic inflammation, plasma cells, and lining hyperplasia. Analysis of RNA sequencing identified 915 significantly differentially expressed genes in correlation with nuclei density (false discovery rate is less than 0.05). Mean nuclei density was significantly higher (P < 0.05) in patients with elevated levels of C-reactive protein, erythrocyte sedimentation rate, rheumatoid factor, and cyclized citrullinated protein antibody. CONCLUSION: Nuclei density is a robust measurement of inflammatory burden in RA and correlates with multiple orthogonal measurements of inflammation.

DRUL for school: Opening Pre-K with safe, simple, sensitive saliva testing for SARS-CoV-2.

To address the need for simple, safe, sensitive, and scalable SARS-CoV-2 tests, we validated and implemented a PCR test that uses a saliva collection kit use at home. Individuals self-collected 300 µl saliva in vials containing Darnell Rockefeller University Laboratory (DRUL) buffer and extracted RNA was assayed by RT-PCR (the DRUL saliva assay). The limit of detection was confirmed to be 1 viral copy/µl in 20 of 20 replicate extractions. Viral RNA was stable in DRUL buffer at room temperature up to seven days after sample collection, and safety studies demonstrated that DRUL buffer immediately inactivated virus at concentrations up to 2.75x106 PFU/ml. Results from SARS-CoV-2 positive nasopharyngeal (NP) swab samples collected in viral transport media and assayed with a standard FDA Emergency Use Authorization (EUA) test were highly correlated with samples placed in DRUL buffer. Direct comparison of results from 162 individuals tested by FDA EUA oropharyngeal (OP) or NP swabs with co-collected saliva samples identified four otherwise unidentified positive cases in DRUL buffer. Over six months, we collected 3,724 samples from individuals ranging from 3 months to 92 years of age. This included collecting weekly samples over 10 weeks from teachers, children, and parents from a pre-school program, which allowed its safe reopening while at-risk pods were quarantined. In sum, we validated a simple, sensitive, stable, and safe PCR-based test using a self-collected saliva sample as a valuable tool for clinical diagnosis and screening at workplaces and schools.

Machine learning integration of scleroderma histology and gene expression identifies fibroblast polarisation as a hallmark of clinical severity and improvement.

OBJECTIVE: We sought to determine histologic and gene expression features of clinical improvement in early diffuse cutaneous systemic sclerosis (dcSSc; scleroderma). METHODS: Fifty-eight forearm biopsies were evaluated from 26 individuals with dcSSc in two clinical trials. Histologic/immunophenotypic assessments of global severity, alpha-smooth muscle actin (aSMA), CD34, collagen, inflammatory infiltrate, follicles and thickness were compared with gene expression and clinical data. Support vector machine learning was performed using scleroderma gene expression subset (normal-like, fibroproliferative, inflammatory) as classifiers and histology scores as inputs. Comparison of w-vector mean absolute weights was used to identify histologic features most predictive of gene expression subset. We then tested for differential gene expression according to histologic severity and compared those with clinical improvement (according to the Combined Response Index in Systemic Sclerosis). RESULTS: aSMA was highest and CD34 lowest in samples with highest local Modified Rodnan Skin Score. CD34 and aSMA changed significantly from baseline to 52 weeks in clinical improvers. CD34 and aSMA were the strongest predictors of gene expression subset, with highest CD34 staining in the normal-like subset (p<0.001) and highest aSMA staining in the inflammatory subset (p=0.016). Analysis of gene expression according to CD34 and aSMA binarised scores identified a 47-gene fibroblast polarisation signature that decreases over time only in improvers (vs non-improvers). Pathway analysis of these genes identified gene expression signatures of inflammatory fibroblasts. CONCLUSION: CD34 and aSMA stains describe distinct fibroblast polarisation states, are associated with gene expression subsets and clinical assessments, and may be useful biomarkers of clinical severity and improvement in dcSSc.

RNA Identification of PRIME Cells Predicting Rheumatoid Arthritis Flares.

BACKGROUND: Rheumatoid arthritis, like many inflammatory diseases, is characterized by episodes of quiescence and exacerbation (flares). The molecular events leading to flares are unknown. METHODS: We established a clinical and technical protocol for repeated home collection of blood in patients with rheumatoid arthritis to allow for longitudinal RNA sequencing (RNA-seq). Specimens were obtained from 364 time points during eight flares over a period of 4 years in our index patient, as well as from 235 time points during flares in three additional patients. We identified transcripts that were differentially expressed before flares and compared these with data from synovial single-cell RNA-seq. Flow cytometry and sorted-blood-cell RNA-seq in additional patients were used to validate the findings. RESULTS: Consistent changes were observed in blood transcriptional profiles 1 to 2 weeks before a rheumatoid arthritis flare. B-cell activation was followed by expansion of circulating CD45-CD31-PDPN+ preinflammatory mesenchymal, or PRIME, cells in the blood from patients with rheumatoid arthritis; these cells shared features of inflammatory synovial fibroblasts. Levels of circulating PRIME cells decreased during flares in all 4 patients, and flow cytometry and sorted-cell RNA-seq confirmed the presence of PRIME cells in 19 additional patients with rheumatoid arthritis. CONCLUSIONS: Longitudinal genomic analysis of rheumatoid arthritis flares revealed PRIME cells in the blood during the period before a flare and suggested a model in which these cells become activated by B cells in the weeks before a flare and subsequently migrate out of the blood into the synovium. (Funded by the National Institutes of Health and others.).

Tranexamic Acid Does Not Reduce the Risk of Transfusion in Rheumatoid Arthritis Patients Undergoing Total Joint Arthroplasty.

BACKGROUND: Patients with rheumatoid arthritis (RA) receive transfusions more often than patients with osteoarthritis following lower extremity total joint arthroplasty (TJA), but mitigating factors are not described. Tranexamic acid (TXA) is widely used to reduce blood loss in patients undergoing TJA, but its effect on transfusion rates in patients with RA has not been studied. METHODS: We retrospectively reviewed data from a prospectively collected cohort of patients with RA undergoing TJA. Disease activity measured by Clinical Disease Activity Index, patient-reported outcome measures, and serologies was obtained. Baseline characteristics were summarized and compared. Transfusion requirements and TXA usage were obtained from chart review. Logistic regression was used to determine factors associated with transfusion in RA patients undergoing TJA. RESULTS: The cohort included 252 patients, mostly women with longstanding RA and end-stage arthritis requiring TJA. In multivariate analysis, 1 g/dL decrease in baseline hemoglobin (odds ratio [OR] = 0.394, 95% confidence interval [CI] [0.232, 0.669], P = .001), 1-minute increase in surgical duration (OR = 1.022, 95% CI [1.008, 1.037], P = .003), and 1-point increase in Clinical Disease Activity Index (OR = 1.079, 95% CI [1.001, 1.162]) were associated with increased risk of transfusion. TXA use was not associated with decreased risk of postoperative transfusion. CONCLUSIONS: Preoperative health optimization should include assessment and treatment of anemia in RA patients before TJA, as preoperative hemoglobin level is the main risk factor for postoperative transfusion. Increased disease activity and increased surgical time were independent risk factors for postoperative transfusion but are less modifiable. While TXA did not decrease transfusion risk in this population, a prospective trial is needed to confirm this. LEVEL OF EVIDENCE: IV.

Rheumatoid Arthritis Flares After Total Hip and Total Knee Arthroplasty: Outcomes at One Year.

OBJECTIVE: Most patients with rheumatoid arthritis (RA) undergoing total hip arthroplasty (THA) and total knee arthroplasty (TKA) have active RA and report postoperative flares; whether RA disease activity or flares increase the risk of worse pain and function scores 1 year later is unknown. METHODS: Patients with RA were enrolled before THA/TKA. Patient-reported outcomes, including the Hip disability and Osteoarthritis Outcome Score (HOOS)/Knee Injury and Osteoarthritis Outcome Score (KOOS) and physician assessments of disease characteristics and activity (Disease Activity Score in 28 joints [DAS28] and Clinical Disease Activity Index), were collected before surgery. Patient-reported outcomes were repeated at 1 year. Postoperative flares were identified using the RA Flare Questionnaire weekly for 6 weeks and were defined by concordance between patient report plus physician assessment. We compared baseline characteristics and HOOS/KOOS scores using 2-sample t-test/Wilcoxon's rank sum test as well as chi-square/Fisher's exact tests. We used multivariate linear and logistic regression to determine the association of baseline characteristics, disease activity, and flares with 1-year outcomes. RESULTS: One-year HOOS/KOOS scores were available for 122 patients (56 with THA and 66 with TKA). Although HOOS/KOOS pain was worse for patients who experienced a flare within 6 weeks of surgery, absolute improvement was not different. In multivariable models, baseline DAS28 predicted 1-year HOOS/KOOS pain and function; each 1-unit increase in DAS28 worsened 1-year pain by 2.41 (SE 1.05; P = 0.02) and 1-year function by 4.96 (SE 1.17; P = 0.0001). Postoperative flares were not independent risk factors for pain or function scores. CONCLUSION: Higher disease activity increased the risk of worse pain and function 1 year after arthroplasty, but postoperative flares did not.